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Using a proposed numerical technique for calculating anomalous Green's functions characteristic 
of superconductivity, we show that the low-lying excitations in a wide parameter and doping region of 
the two-dimensional t— J model are well described by the picture of dressed Bogoliubov quasiparticles 
in the BCS pairing theory. The pairing occurs predominantly in d x 2_ y 2-wa,ve channel and the energy 
gap has a size Ad— 0.15 J— 0.27 J between quarter and half fillings. Opening of the superconducting 
gap in the photoemission and inverse-photoemission spectrum is demonstrated. 
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Recent issues in the theory of high-temperature super- 
conductivity include the fundamental question whether 
the two-dimensional (2D) t—J model contains a super- 
conducting phase relevant to cuprate materials, and if so, 
what types of superconductivity the model has. There 
are a number of mean-field based, variational Monte 
Carlo, and other approximate calculations, but a com- 
monly agreeable phase diagram has still been out of our 
reach. Thus, at this stage, unbiased numerically-exact 
calculations of the model certainly offer important infor- 
mation on this question. Among such exact calculations 
, Dagotto and Riera have recently found indica- 
tions of superconductivity near the region of phase sep- 
aration (J] by examining equal-time pairing correlations 
in finite-size clusters; the existence of a ci^-j^-wave con- 
densate has thereby been conjectured. 

In this Letter, we propose a new technique for examin- 
ing the low- lying excitation spectrum, i.e., the exact cal- 
culation of anomalous Green's functions for Bogoliubov 
quasiparticles in finite-size clusters; we can thereby calcu- 
late the quasiparticle excitation spectrum directly to ex- 
amine the superconducting pairing interactions in the 2D 
t—J model. Then, we show that the low-energy excita- 
tions in a wide parameter ( J/t) and doping region of the 
model are well described by the picture of dressed Bogoli- 
ubov quasiparticles within the BCS (Bardeen-Cooper- 
Schrieffer) pairing theory: the singlet pairing of electrons 
with opposite momenta occurs near the Fermi energy un- 
like in the regime of Bose condensation of a gas of spin 
singlets in real space. The pairing occurs predominantly 
in the c? a ,2_ l) 2-wave channel, and the gap energy scales 
with J and has a magnitude Ad~0.15 J— 0.27J. We also 
demonstrate the opening of a superconducting gap in 
the photoemission and inverse-photoemission spectrum. 
This work provides the first demonstration of the valid- 
ity of the BCS pairing theory for a strongly-correlated 
electron model. 

The t—J model is defined by the Hamiltonian 
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where c\ a (ci a ) is the projected electron-creation (annhi- 
lation) operator at site i and spin a (=T>I) allowing no 
doubly occupied sites, Si is the electronic spin operator, 
and rii= njf+n^ is the electron number operator. The 
summation is taken over all the nearest-neighbor pairs 
<ij> on the two-dimensional square lattice. The aver- 
age number of electrons per site, n, is defined with the 
terms half filling (n=l) and quarter filling (n=l/2). 

We adopt the working hypothesis that low-lying states 
of the clusters can be described by the microcanonical 
version of the BCS pairing theory ||. Hence, we assume 
that a cluster ground-state with an even electron-number 
N can be written as \i(}^)=P N \BCS) where |BCS) is a 
BCS wave function (formulated in terms of 'quasipar- 
ticles' so that the effects of strong correlations are al- 
ready included) and P/v projects on the TV-electron sub- 
space. Similarly, we assume that the low-lying states 
with an odd electron- number N+l can be written as 
|-0^ r+1 )=i 3 jV"-l-i'y^ a .|BCS) where 7 kcr creates a Bogoliubov 
quasiparticle with momentum k. Then, if we define the 
following anomalous Green's functions, our working hy- 
pothesis predicts that they should describe the excita- 
tions of Bogoliubov quasiparticles in the clusters. 

The one-particle anomalous Green's function is defined 

as 

G(k,z) = « +2 l4 T z _ g 1 +£ , ^-kjO 

where c kcr is the Fourier transform of cj CT and Eq is the 
ground-state energy. We define the spectral function 
F(k,u;)=— (l/TrjlmG^k, u + irj) with rj—0 + and its fre- 
quency integral -Fic=(V'o' +2 | c kT c -k].l' ! /'o r )- The hypothesis 
then predicts F(k, u>)=Fi i 5(uj — Ek) with F k =z k A k /2i?k: 
via the Bogoliubov- Valatin transformation of operators 
[^J, where Zk, i?k, and Ak are the wave- function renor- 
malization constant, renormalized quasiparticle energy, 
and gap function, respectively. Thus, this Green's func- 
tion describes the excitation of one Bogoliubov quasipar- 
ticle in the cluster. Similarly, we define the two-particle 
anomalous Green's function as 



G(k,k',z) = (V>o r |cl k , ; c_ k .i — =]y C k , T Cktl-0o 
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with rikcr=(V'o'l c kcr c kcr|V'o r ) an d the ground-state energy 
Eq . Only the TV-particle subspace is involved unlike 
in Eq. (Q). We define the spectral function F(k, k',w) 
and its frequency integral -fkk' as above. The hypoth- 
esis then predicts F(\c, k', ui)=F^F^,5(uj~E] il ~ 2?k') : the 
Green's function describes the excitation of two Bogoli- 
ubov quasiparticles. Note that Eq. (||) may also be de- 
fined with c k , CT and Ck<j interchanged and with nko- re- 
placed by nko-=(V'o r | c ko-c ko .| , 0j^); whereas no change oc- 
curs in the BCS theory, some difference appears in the 
calculated spectra because of the violation of the anti- 
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commutation relation of the constrained operators; how- 
ever, we find the difference to be insignificant. Thus, 
by examining these two anomalous Green's functions, we 
can see if the low-energy excitations of the 2D t— J model 
are described by the BCS pairing theory. 

The off-diagonal Green's functions Eqs. (0) and (j|) 
may be evaluated by subtraction of diagonal ones, which 
in turn are computed in the standard Lanczos algorithm 
0: e.g., for the one-particle anomalous Green's function, 

we prepare the state CktlV'(D /+2 )+ c -kjjV'o V ) anci calculate 
the spectral function, from which the usual single-particle 
spectral functions for Ckt IV'o^ 2 ) and c_ k ,\ipQ) are sub- 
tracted; thereby F(\c,u>)+c.c. is obtained. The majority 
of the spectral weight is cancelled out, remaining only the 
anomalous part (compare the spectra shown below) . The 
imaginary part of F(k, to) is also calculated by using the 
state Ck||^/'(^ r+2 )+icL k J'(/'o r ) an d is always found to van- 
ish. We use clusters of the size 4x4 and VT8x vl8 with 
periodic boundary condition; the Green's functions are 
evaluated for zero-momentum ground states to simulate 
the excitations in the thermodynamic limit. The Fermi 
momentum kp is located at the fc-points where F(k, w) 
and F(k, k', uS) show the lowest-energy peak. 

The calculated results for F(k, w) are shown in Fig. [j]. 
We choose the average value of Eg and E^ +2 as the 
value of -Bo used in Eq. (|j). We find the following, all 
of which are consistent with expectations of the BCS 
pairing theory: A pronounced low-energy peak appears 
at kp and smaller peaks appear at higher energies for 
other momenta; the weights of the peaks are consistent 
with the BCS form of the condensation amplitude Fk 
with a maximum at kp. The momentum dependence of 
F(k, w), i.e., the change in sign under rotation by tt/2 
and vanishing weight along the k x =k y line, is a clear in- 
dication of ci^ _ y 2 -wave pairing. The size of the energy 
gap, which may be estimated directly from the positions 
of the peaks, scales well with J/t value. With decreasing 
gap size, the peaks at momenta other than kp lose their 
weight as expected from the BCS theory. An important 
consequence of d x 2_ y 2-wave pairing may be seen in the 
point-group symmetry of the ground states, i.e., an alter- 
nation of the symmetry between A\ and B\ in, e.g., the 
18-site cluster for fillings of 10, 12, 14, and 16 electrons 
in a fairly wide J/t region. This alternation, absent in, 
e.g., the attractive-C/ Hubbard clusters, suggests the pic- 
ture that electrons are added in pairs with d :c 2_j / 2-wave 
symmetry. G(k, z) in Eq. ([[]) (and thus F(k, u>)) reflects 
the pairing symmetry via the point-group symmetries of 
■>Pq and '0o V+2 - Note that F(k, k',w) is defined entirely 
within the A^-particle subspace and thus is not affected 
by this alternation; nevertheless it shows the same indi- 
cation of <i 2 .2_ y 2-wave pairing as F(k, uj). 

The calculated results for F(k, k',o;) are shown in 
Figs. [2] and ||. We again find that the size of the excita- 
tion gap scales well with J/t, and the k dependence of 
the spectra clearly indicates d x 2_ y 2--wave pairing. There 
are sharp peaks at low energies and broadened features 
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at higher energies. As is the case for F(k, to), these high- 
energy features lose their weight rapidly with decreasing 
J/t value, whereas the peaks at kp become sharp but re- 
main finite with decreasing J/t value. These results are 
consistent with the notion of weakly-interacting Bogoli- 
ubov quasiparticles for low-lying excitations in the BCS 
superconductors. 

To be more quantitative, let us compare the cal- 
culated spectra F(k, w) and F(k, k',w) with the BCS 
predictions for the spectral functions: the quasipar- 
ticle energy £ic=\/£k + A£ with one-electron energy 
£k=— 2t e g (cos k x + cosfcj,)— ji and gap function Ak = 
A,j(cos k x — cos k y ) is assumed. We use the effective hop- 
ping parameter t c g to take into account the quasiparti- 
cle band narrowing. The chemical potential /i is chosen 
to guarantee vanishing £k at kp. The value of A<j is 
then evaluated by fitting the positions of the low-energy 
peaks in -F(k, to) and F(k, k', u>). The renormalization 
factor zb (defined by assuming z^ to be k-independent) 
is estimated from weights of the low-energy peaks in 
F(k,uj) and -F(k, k',w). The parameters t e g and z-q re- 
flects the effects of strong correlation and imply the use 
of a dressed Bogoliubov-quasiparticle description. The 
fitted quasiparticle spectra are shown by dotted curves 
in Figs. [j]-^ We find a fair agreement with the ex- 
act spectra, which demonstrates the validity of the BCS 
pairing theory for low- lying excitations in the 2D t— J 
model. This is also the case at low-doping levels: both 
F(k,uj) and F(k,k',w) exhibit well-defined low-energy 
peaks with d 2 .2_ a 2-wave symmetry which can be fitted to 
the BCS spectra with rather small t e g and zb values. 

The estimated values of the gap parameter A^ are 
shown in Fig. [|. We find that A^ scales well with 
J and has the value Ad~0.15J— 0.27J until reaching 
the region of phase separation. The gap value is large 
(<0.27J) near half filling and small (~0.15J) around 
quarter filling. The renormalization factor varies over 
Zb— 0.2— 1 with smaller values at low-doping levels. Note 
that the maximum gap-energy 2Arf/i~0.9 (0.5) around 
quarter filling at J/t— 2.5 (1.5) (see Fig. ||) is significantly 
smaller than the effective half-bandwidth 4t c g /t~2.2 esti- 
mated from the fitting of F(k, u) and F(k, k',w). Also, 
at low-doping levels, we note that the estimated value 
2Ad— 0.5 J— 0.6 J is smaller than the effective bandwidth 
of ~2 J— 4J. Thus, even near phase separation or at low- 
doping levels, the superconductivity in the 2D t— J model 
is not in the regime of Bose condensation of a gas of spin 
singlets in real space. 

The single-particle spectral function A(k, u>) (see Ref. 
for the definition) simulates angle-resolved photoe- 
mission (PES) and inverse photoemission (IPES) spec- 
troscopy. The calculated spectra are shown in Fig. ||, 
which shows how the superconductivity affects the single- 
particle excitations. We find that with increasing J/t a 
gap-like feature appears at k=(27r/3,0) which may be 
interpreted as a superconducting gap as the BCS disper- 
sion ±_Ek suggests. An associated spectral-weight trans- 
fer also indicates a tendency toward smearing of the jump 
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at kp in the momentum distribution function. We have 
also calculated the spin-excitation spectrum S(q, oj) (not 
shown) and compared it with the bare susceptibility in 
the BCS mean-field theory. We find overall agreement in- 
cluding spectral-weight distributions and opening of the 
superconducting gap: the features are well interpreted in 
terms of the particle-hole excitations across the nonin- 
teracting Fermi surface although at q=(7r,7r) the effect 
of nearest-neighbor antiferromagnetic spin correlations 
becomes appreciable in the large- J/t region or at low- 
doping levels. Details will be discussed elsewhere. 

Summarizing, we have studied the superconductivity 
in the 2D t—J model with the use of a newly proposed 
technique for examining the low-lying excitation spec- 
trum, i.e., the exact calculation of anomalous Green's 
functions for Bogoliubov-quasiparticle excitations char- 
acteristic of the superconducting state. The validity of 
the BCS pairing theory for a strongly-correlated electron 
model has thereby been demonstrated for the first time. 
We have shown that the pairing occurs predominantly in 
the d x 2_ y 2 -wave channel and the energy gap has a size 
A d ~0.15J— 0.27 J in a wide region between quarter and 
half fillings. 
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FIG. 1. Bogoliubov-quasiparticle spectrum F(k,u>) for 
the 16-site cluster with doping between 4 and 6 holes at (a) 
J/t=0.8 and (b) 1.5, and for the 18-site cluster with dop- 
ing between 4 and 6 holes at (c) J/t=0A and (d) 1.5. The 
spectra at other momenta k are identical with those shown 
here except the sign that follows the d x 2_ y 2-wave symmetry; 
the spectra along the line k x —k y vanish exactly. kF is lo- 
cated at (-7T, 0) in (a) and (b), and at (27r/3, 0) in (c) and (d). 
The dotted curves show the BCS spectral function obtained 
for parameter values Ad/t=0.13 in (a) and 0.29 in (b) with 
t eff /i=0.55 and z B = 1.0, and A d /£=0.05 in (c) and 0.27 in 
(d) with t e ff/t=0.45 and zb=0.5. We use the value r?/t=0.15. 



FIG. 2. Bogoliubov-quasiparticle spectrum F(k, k',a>) for 
the 16-site cluster with 8 holes at J/t=0A (left panel) and 
2.5 (right panel). The momentum k dependence is shown 
with k' fixed at (0, — tt/2). kF is located at (7r/2,0) and its 
equivalent points. The value »?/i=0.15 is used. The dotted 
curves show the BCS spectral function obtained for parameter 
values A,j/t=0.105 (left panel) and 0.47 (right panel) with 
t e g/t=0.55 and zb=0.7. 

FIG. 3. As in Fig. ||but for the 18-site cluster with 6 holes 
at J/t=0.4 (left panel) and 1.5 (right panel). The momentum 
k' is fixed at (0, — 27r/3), and kF is located at (2w/3, 0) and 
its equivalent points. The value r//t=0.15 is used. The BCS 
spectra are for the parameter values Ad/£=0.05 (left panel) 
and 0.285 (right panel) with t e s/t=0A5 and zb=0.5. 

FIG. 4. Gap parameter Ad for various J/t values and 
doping rates. The panel gives the number of holes (in the 16- 
or 18-site cluster) in the final state of the anomalous Green's 
functions: odd (even) numbers imply that the results are from 
the one-particle (two-particle) anomalous Green's function. 
The dashed lines are a guide to the eye. 

FIG. 5. Single-particle spectral function A(k, uS) for the 
18-site cluster with 6 holes at J/t=0A (left panel) and 1.5 
(right panel). u> is measured from the 6-hole ground-state 
energy, k's are shown in the panels, and the value »7=0.05 
is used. The thin-dashed curves show the BCS dispersion 
±-Bk obtained for the parameter values used in Fig. Bp good 
agreement is seen near the Fermi energy while at higher en- 
ergies meaningful comparison cannot be made because of the 
absence of damping effects. 
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